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Abstract. We establish unique existence and stability of subsonic potential flow for steady Euler- 
Poisson system in a multidimensional nozzle of a finite length when prescribing the electric potential 
difference on non-insulated boundary from a fixed point at the exit, and prescribing the pressure 
at the exit of the nozzle. The Euler-Poisson system for subsonic potential flow can be reduced to a 
q , nonlinear elliptic system of second order. In this paper, we develop a technique to achieve a priori 

C 1 ' 01 estimates of solutions to a quasi-linear second order elliptic system with mixed boundary 
conditions in a multidimensional domain with Lipschitz continuous boundary. Particularly, we 
discovered a special structure of the Euler-Poisson system which enables us to obtain C 1 '" estimates 
of velocity potential and electric potential functions, and this leads us to establish structural stability 
of subsonic flows for the Euler-Poisson system under perturbations of various data. 

< 

^ \ 1. Introduction and Main Results 

1.1. Preliminary. The following nonlinear system, called the Euler-Poisson system, 

' p t + div(pu) = 0, 
(pu) t + div(pu ® u + pin) = pV$, 
(p£)t + div(pf u + pu) = pu ■ V$, 
A$ = p - b(x), 



(N 



>: (i.i) 

m 

models various physical phenomena including the propagation of electrons in submicron semicon- 
ductor devices and plasmas (cf.|22j). and the biological transport of ions for channel proteins (cf.[2]). 
In the hydro dynamical model of semiconductor devices or plasmas, u, p, p, and £ represent the 
macroscopic particle velocity, electron density, pressure, and the total energy, respectively. The 
electric potential $ is generated by the Coulomb force of particles. I n is the n x n identity matrix 
and b(x) > stands for the density of fixed, positively charged background ions. In the biological 
model describing transport of ions between the extracellular side and the cytoplasmic side of the 
membranes(cf. [2]), p, pu and $ are the ion concentration, the ion translational mass, and the 
electric potential, respectively. The system (jl.ip is closed with the aid of definition of total energy 
and the equation of state 

lul 2 

(1.2) £ = U- + e , p = p(p,e), 

where e is the internal energy. By the second law of thermodynamics, the entropy S is a constant 
along each particle trajectory in smooth steady flow of (jl.l|) . So the entropy S is globally a constant 
if S is a constant at the entrance of flow. In this paper, we consider the case when the entropy S 
is globally a constant, and such a case is called isentropic. For isentropic flow, we can write the 
pressure p and the internal energy e as 

n o, f p p'(e), p(p) 

(1.3) P=P{P), e= / dg 

Jko 6 P 
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for a constant k$ > 0. We assume that p € C°([0, +00)) n C 3 ((0, 00)) satisfies: 

(1.4) p(0) = 0, > 0, p"(p) > for all p > and p(+oo) = +00. 



For a constant 7 > 1, p(p) = p' 1 is a typical example of p satisfying f 1 1 . 4 [) . which corresponds to 
polytropic gas in gas dynamics. 

The steady Euler-Poisson system of isentropic flow is 



for a given function b > 0. There are several issues to make the system (|1.5p complicated. The 
first is that fll .5f) is a mixed type system, and its type depends on the Mach number M which is 
given by M = ^jy for c(p) = \fp'{p). Here, c is called the local sound speed. If M < 1, then (|1.5p 
can be decomposed into a nonlinear elliptic system and homogeneous transport equations, and the 
flow is said subsonic. If M > 1, on the other hand, (II. 5j) can be decomposed into a nonlinear 
hyperbolic-elliptic coupled system and homogeneous transport equations at best, and the flow is 
said supersonic. The second issue is that the last equation in fll .5f) . which is a Poisson equation, has 
a nonlocal effect to the other equations in (II. 5ft . and it makes the fluid variables p, u and electric 
potential $ interact in a highly nonlinear way. Also, physical boundary conditions such as fixed 
exit pressure give nonlinear boundary conditions for the system (|l,5p . 

It is our goal to prove unique existence and stability of subsonic flows for steady Euler-Poisson 
system in a multidimensional nozzle under perturbations of exit pressure, electric potential differ- 
ence on non-insulated boundary and under perturbations of the nozzle itself. Our motivation is 
originated from the study on structural stability of transonic shocks occurring in flow governed by 
steady Euler-Poisson system. There have been a few studies on transonic shocks of Euler-Poisson 
svstemfcf . [T| \9[ \19 \ \20 [ [25]) . In [20], the authors considered one dimensional solutions of (jl.5p with 
a constant background charge b(x) = 60 > 0, and proved the unique existence of transonic shock 
solutions provided that the entrance and exit data are appropriately given. So it is natural to 
ask whether these one dimensional transonic shocks are dynamically and structurally stable. The 
dynamical stability of one dimensional transonic shock solution is achieved in [19] . In order to 
study structural stability of transonic shocks, however, it is inevitable to consider small pertur- 
bations of one dimensional transonic shocks in multidimensional domain, but there are very few 
known results about multidimensional solutions of Euler-Poisson system(cf.[9l [TO]). Comparing 
with extensive studies and recent significant progress on transonic shock solutions of the Euler sys- 
tem(see [31 [6j [26] and references therein), stability problems for multidimensional transonic flows 
of the Euler-Poisson system are essentially open. The main difference of the Euler-Poisson system 
from the Euler system is that the Poisson equation for electric potential is coupled with the other 
equations in the Euler-Poisson system. And, this makes it hard to analyze even one dimensional 
solution of the Euler-Poisson system. In fact, one dimensional flow of the Euler-Poisson system 
behaves very differently from the one of the Euler system(cf.[20j). And it is even harder to study 
multidimensional transonic flow of the Euler-Poisson system due to nonlinear interaction between 
the electric potential and all the other fluid variables. 

As the first step to investigate stability of multidimensional transonic flow of the Euler-Poisson 
system, we establish the unique existence and stability of subsonic flows of steady Euler-Poisson 
system under perturbations of the exit pressure and electric potential difference on non-insulated 
boundary. There have been a few results about existence of subsonic solution of hydrodynamic 
equations, which are the Euler-Poisson system with relaxation, under smallness assumptions on 
the flow velocity for both unsteady and steady cases(see El [121 [21]). In this paper, we prove 
existence of multidimensional subsonic solutions of steady Euler-Poisson system without smallness 
assumption on the flow velocity. 



div(pu) = 0, 



(1.5) 



< 




A$ = p - b 
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Fix an open, bounded and connected set A C M" _1 (n > 2) with a smooth boundary dA, and 
define a nozzle M by 

(1.6) M := {x = (x, x n ) € R n : x' e A, x n E (0, L)} C M n . 

The nozzle boundary dJ\f consists of the entrance To = A x {0}, the exit Tl = A x {L}, and the 
insulated boundary T w = dA x (0,L). In order to study the system (|1.5fl in a multidimensional 
domains TV with arbitrary cross-section A, we consider irrotational flow where the velocity u of the 
flow is represented by 

(1.7) u = Vp 

for a scalar function cp which is called a velocity potential function. By (11.31) and (jl.7p . the second 
equation in (11.5P can be rewritten as 

(1.8) /A7(^-$) = 
for 

(1-9) ^=:W + r 

For p > 0, (fl~8j) implies 

(1.10) &-$ = K 

for a constant i^o- Without loss of generality, we choose Kq = 0. Set 

(1.11) h{p) := r^-dg. 



Jk a 

Then, the equation (jl.lOp with Kq = implies that h(p) = $ — ^|V<p| 2 . From the condition (jl.4p . 
we have h'(p) > for all p > and lim hip) = oo. So the function h is invertible wherever 

p— s-oo 

$ — ^|V<y?| 2 > liminf h(p), in which case one can rewrite (jl.lip as p = /i _1 (<I> — ^|V(^| 2 ). We use 
this expression to reduce (|1.5|) to a nonlinear system of second order equations for <p and $: 

(1.12) div(p($, |Vv?| 2 )V(p) = 0, 

(1.13) A$ = p($, |V99[ 2 ) - b 
with p > given by 

(1.14) p{$>,\V^) = h-\$>- l -\Vv?) 

provided that h^ 1 is well defined. If we regard (|1.12p as an equation for ip, then it is mixed type. 
More precisely, (|1.12j) is elliptic if and only if 

(1.15) |V<p| 2 < p'(p) (subsonic) 
and is hyperbolic if and only if 

(1.16) \V <p\ 2 > p' (p) (supersonic) . 

The system of (|1.12p and (|1.13p becomes a quasilinear elliptic system if (|1.15p holds, and a 
hyperbolic-elliptic coupled system if (|1.16p holds. 

Our interest is on stability of subsonic solution under perturbations of exit pressure and electric 
potential difference on non-insulated boundary from a fixed point. So the boundary conditions are 
formulated as follows. First, for a given function p ex on Tl, set 

(1.17) p(p(<S>,\V^\ 2 )) =p ex on T L . 

On the wall boundary, slip/insulated boundary conditions for ip and are prescribed as follows: 

(Lis) d nwV = d nvj $ = o on r„ 
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where n„, is the unit inward normal vector on r^. We fix a point xq on T^, and prescribe the 
electric potential difference between two points x G Tq U Yl and xq as follows: 



(1.19) $(x) - $(xo) 



$ en (x) for x G T 
$>ex(x) for x G T L 



In (I1.19p . the value of 3>(xo) is uniquely determined by (jl.lOp and one point boundary condition 
for the Bernoulli's function: 

(1.20) #(xo) = ^o- 

Remark 1.1. Differently from the Euler system, the Bernoulli's function 3$ = §|VV?| 2 + /^ ^—^-dg 
is not a constant along each streamline due to the equation (|1 .8j) . So we call SS the Bernoulli's 
function rather than the Bernoulli's invariant. This is one of differences of the Euler-P oisson 
system from the Euler system. 

Finally, homogeneous Dirichlet boundary condition for tp is imposed at the entrance: 

(1.21) ip = on T . 

The boundary data p ex , $ en , <& ex and 9S§ will be specified later. 

1.2. Main theorem. If we fix b as a constant 6q > in the equation ()1.13j) then the equations 
(|1.12p and (11 . 13f) become invariant under translation. So if the boundary data $ en , Q ex and p ex 
are all constants, then one may look for a solution (y>, $) as functions of x n only for x n G (0,-L). 
We note that $ ex = if $ ex is a constant. 

Proposition 1.2 (One dimensional subsonic flow). Fix constants bo > and L > 0. TTien 
i/iere exists a nonempty set *Po o/ parameters in R 2 x R+ so i/iat /or any ($ e n,0) ^0,0? Pes, o) G 
*Po) */ (^enj ^dPea;) = (^enfi, 33o,o, Pex,o) then the system of (I1.12P and (I1.13P in A/" with the 
boundary conditions 17f l f)1.21 j) has a unique C 2 one- dimensional solution (cp, $) inJ7 satisfying 
the inequalities p($, \V(p\ 2 ) > and |Vi^| 2 < p'(p($, |V^| 2 )) in AT. 

The proof of Proposition 11.21 is in Section 11.31 

Remark 1.3. We can find one- dimensional solutions for a nonconstant function b provided that 
the boundary data (J&en,0i ^o,o>Pea;,o) are properly chosen depending on b. Details can be found in 
Appendix 0J 

Definition 1.4 (Background solution). Fix constants bo > and L > 0, and let the parameter 
set ^Po be as in Proposition For given constants (& e n,o, &o,o,Pex,o) £ tyo, let &)(%n) be the 
corresponding one- dimensional solution to the system of (|1.12p and (|1.13p in N with the boundary 
conditions f| 1 . X T j) — (| X . 2 1 j) for the boundary data ($ e n,o> ^o,o, Pex,o)- For such solution (cp, <J>)(x n ), we 
define two functions (ipo, &q)(x\, • • • , x n ) in N by 

<Po(x lr -- ,x n ) := <p(x n ), 
$ Q (xi,--- ,x n ) := $(x n ) 



We call (ipo,$>o) the background solution corresponding to (& e n,o, 3$o,OiPex,o) £ ^Po- and &o ire 
smooth in U and satisfy 

(1.23) p($o,|V^ | 2 ) > ^o, p'(p($o,|V^o| 2 ))-|V99 | 2 >^o in AT 

for a constant vq > depending on bo, L, & ex , £§o,o an d Pex,o- 

Our goal is to prove stability of a background solution under small perturbations of the back- 
ground charge b and the boundary data. In Appendix [HI we also consider stability of a background 
solution under small perturbations of the nozzle J\f. 
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Problem 1. Fix bo > and L > 0, and let the parameter set *}3o be as in Proposition 
Given ($ e n,o> ^o,o,Pex,o) £ *Po; (VOj3?o) be the corresponding background solution. Also, let 
($ en , &ex,Pex) be small perturbations of ($ e n,o 5 ®,Pex,o) and let a constant £$o be close to &o,o- 
And, let p(&, |V(/?| 2 ) > be defined by (|1.14p , Find a solution ((/?, <I>) to i/te nonlinear boundary 
value problem: 

(1.24) |div(^|V r P)V,) = 

(1.25) <^ = 0onr , p(p(<S>,\V<p\ 2 )) = Pex on F L 

(1.26) * = + ° nr ° 

{<%o + $ ex on T L 

(1.27) d n ^ = d nm $ = on r tt . 

In order to state our main results on Problem [TJ weighted Holder norms are introduced first. 
For a bounded connected open set fi C R n , let T be a closed portion of d£l. For x,y£S], set 

<5 X := dist(x, T) and 5 X)J := min(<5 x , <5 y ). 

For k G R, a G (0, 1) and m G Z + , define the standard Holder norms by 

„ , , m r i |£>^(x) -Z^u(y)l 
IMIm,o,fi : = 2^ sup | it (x) |, [wJ m ,Q !l n := 2^ SU P i _ i Q > 

0<|p|<m |p|=m 

and the weighted Holder norms by 

HiSi- £ S u P c ax(l/3|+M) i^(x)i, 

0<|/3|<m xefi 

r„l(fe.r) ._ V- rmax(m+a+fc,0) |£>M*) - D & u{j)\ 

m m ,a,n — Z_> S ";P Iy- v I« 

l^l^x.yen.x^y |x y| 

ii n ii n i r i ii n( fc T) n n( fc ' r ) , r i( fc > r ) 

where denotes dt\ ■ ■ ■ d& for a multi-index = (ft, • • • M wi th Pj G Z+ and \/3\ = YTj=x Pj> 
C^'ps(fi) denotes the completion of the set of all smooth functions whose || • ||^'^ q norms are finite. 
The main result of this work is the following. 

Theorem 1. Let J\f be as in (|l,6p . Fix bo > and L > 0, and let the parameter set be as in 
Proposition I l.SX Given (& e n,o, 3$o,o,Pex,o) £ ^Po; tei (^o^o) be the corresponding background solu- 
tion. Assume that b, (& en ,$ ex ,p ex ) and ^o are given as small perturbations of bo, ($ e n,o> ®:Pex,o) 
and ^o,o> respectively, in the following sense: 

\\b - b \\ a ,M < a, 

(1-28) ||$ en - 5>en,o||2,a,ro + II ®ex h,a,T L + \\Pex ~ Pex,o\\a,V L < 

\#0-&Q,0\ <° 

for a small constant a > to be specified below. Also, suppose that <& en and <& ex satisfy the 
compatibility conditions 

(1.29) d nw $ en = on To n T w , 8 nw $ ex = on T L n IV 

Then, /or any <yroen a G (0,1), £/iere exists a constant a > depending on bo, L,& en fi,£$o,OiPex,o 
and a such that wherever a G (0, a], if the boundary data and b satisfy (|1.28|) and (|1.29|) . then the 
nonlinear boundary value problem (|1.24p ~ (|1.27p has a unique solution (ip, $) G [C 1,a (Af)nC 2,a (J\f)] 2 
satisfying the following properties: 
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(a) The equations in (jl.24p form a uniformly elliptic system in M . Equivalently, the solution 
(p, $) satisfies the inequality 

p'(p(<$>, \Vip\ 2 )) - \Vip\ 2 > v > in J7 

for a positive constant v, i.e., the flow governed by (ip,$) is subsonic; 

(b) (p, <3?) satisfy the estimate 

h m\ \\,n m nC-i-^-rouri) , , (-i-^rouri) „ 

for a in (|1.28p . T/ze constants v and C depend only on bo, L,^ en ^, ^Q t Q,p ex fi,n, A and a. 

For simplicity, we hereafter say that a constant C depends on the data if C depends only on 
bo,L,$ enfi ,& 0fi ,pex,o,n,A and a. 

We point out that the boundary conditions (|1.17p - (|1.2ip are physically measurable. In one- 
dimensional solutions, they correspond to prescribing the pressure (or equivalently prescribing the 
density) at both ends of the nozzle. See Remark 11.51 for details. 

In order for a solution (ip, $) of (|1.12p and (11.13P to satisfy fjl . 15|) . it is essential to establish 
L°° estimates of (Dip, $). It is the new feature of this work that we prove C 1,a regularity of 
solutions to a class of second order elliptic systems with mixed boundary conditions on a Lipschitz 
continuous boundary, and use the result to find a solution (ip, <I>) of (|1.12p and (|1.13p so that (ip, <I>) 
satisfy (I1.15P in the nozzle M. Furthermore, the C 1,a estimates that we achieve applies up to the 
boundary dJ\f. Differently from elliptic equations, C 1,a regularity of second order elliptic system is 
not generally known. In the Euler-Poisson system, however, we discovered a special structure(see 
property (c) of Lemma |2 . 2[) . and use it and the divergence structure of the Euler-Poisson system 
to get a priori H 1 estimate of weak solutions to second order elliptic system which yields a priori 
C 1,a estimates of weak solutions to second order elliptic systems with mixed boundary conditions 
in the Lipschitz domain M. Thanks to this property, we can prove unique existence and stability 
of multidimensional subsonic solutions of Euler-Poisson system without smallness assumptions on 
the current flux or flow velocity. This distinguishes the results in this paper from previous works in 
[3 Q21 [21] . Furthermore, the technique developed in this paper can be used to deal with subsonic 
flow of nonisentropoic Euler-Poisson system with nonzero vorticity. This case will be discussed in 
the upcoming papers. 

1.3. One dimensional subsonic flow. Before we proceed further to prove Theorem [lj we first 
prove Proposition 11.21 

Proof of Proposition For simplicity, we regard a one dimensional solution (ip, <3?)(xi, • • • ,x n ) 
in M as functions (ip, &)(x n ) on the interval (0, L). Note that any functions (ip, $)(x n ) satisfy the 
boundary conditions (|1.18p on r^j. Suppose that (ip(x n ),*§*(x n )*) is a one-dimensional solution to 
the system of (|1.12p and f 1 1 . 1 3 1) in M, and let us set u := d Xn ip and E = d Xn <&. Then, for p defined 
by (|1.14p . (p,u,E) satisfy 

(pu)' = 

(1.31) (pu 2 +p(p))' = P E 

E' = p-b 

in the interval (0, L) where ' is the differentiation with respect to x n . The first equation in (11.3ip 
implies that 

(1.32) pu = J on (0,L) 

for a constant Jo > 0. Then by repeating the analysis in [20] for the ODE system (11 .31 H . we can 
state as follows: 



SUBSONIC FLOW FOR MULTIDIMENSIONAL EULER-POISSON SYSTEM 7 

Fix constants 60 > and L > 0. For any given constant Jo > 0, there exists a nonempty set 
*Pi(Jq) C M + x 1 so that for any (po, Eq) G ^Pi(Jq), t/ie initial value problem 



(1.33) 



« = ^ 



(pu 2 +p(p))' = pE on (0,L), 



o 



E' = p-b 

has a unique smooth solution (p,u, E)(x n ) on the interval [0,L] with satisfying 

(1.34) p>0 and p'(p) - u 2 > v x > on [0,L] 

/or some positive constant v\ where the choice of v\ depends on Po,Eq, Jq and L. 

Hereafter, ^i(Jq) denotes the maximal set of (pq,Eq) for which (|1.33p has a unique smooth 
solution (p,u,E) satisfying (11.34p . 

For each Jo > and (po, Eq) E *Pi(Jo), we claim that there exists unique corresponding boundary 
data ($en,0) 3$o,OiPex,o) for boundary conditions (|1.17p - (|1.20p . Fix Jo > and (pq,Eq) G *Pi(Jo). 
The first equation of ()1.33|) yields u(0) = j^. By substituting this into (|1 .9j) . we get ^(0) = 

±(^) 2 + ^fidg. Therefore $(0) = &(0) = \{^f + ^dg is obtained from the equation 
3D — = 0. Using this equation again gives $(£) = =^o,o and <1?(0) = Q e n,o + ^o,o- Therefore, we 
get the following equation for $ en; o and ^o,o 

(1.35) $ enfi + ^o,o =u-) 2 + r — d e- 

2 \PoJ Jko Q 

For the solution (p,E) to (|1.33|) . set 

(1-36) p exfl = p{p{L)). 

Substituting p = p(L) and u{= |V<p|) = into (|1.9p yields 

(i.3T) % . 9(i)= if4V + r^, 



2 \p(L)J 



hi 



By (|1.35p and (11.37p . the value of $ e n.o is uniquely determined by (pq,Eq) and L. For such 
(&en,o, -^o,o,Pex,o), define two functions ip and $ by 

rx„ f-x„ 

(1.38) ip(x n ):= u{t)dt, $(x n ) := ^o,o + ^en.o + / E(t)dt 

Jo Jo 

for the solution (p, .E) to (ll.33p . Then the pair of (cp, 3>) becomes a C 2 solution of (fl~12jl and (fL13l> 
in J\f with the boundary conditions (|1.17p - (|1.2ip . Moreover, by (jl.34p . we have 

(1.39) p'(p($, |V(p| 2 )) - |V<p| 2 > z/i in AT 

for v\ > same as in (ll,34p . This proves the existence of nonempty parameter set *Po- Next, we 
prove the uniqueness of one-dimensional solution for fixed (& e n,o, 3$o,o,Pex,o) £ *Po- 

Fix the boundary data (i erii o, &o,o,Pex,o) G ^Po, and let (<pW, $W) and (<p( 2 ), $( 2 )) be two one- 
dimensional solutions of (|1.12p and (|1.13p in Af with the boundary conditions (|1.17p - (|1.2ip and 
satisfying (fiTT5|) . For each j = 1,2, let be defined by (fl"T4"l) for (tp& , $0')) . By (JOJ), for any 
given p ex> o > 0, there exists unique p e z,o > satisfying p{p ex ,o) = Pez,o- Then each «vi) = (tp(i)y 
satisfies 



(1.40) -{uW{L)Y = &w- ^dg for j = 1,2. 

z Jk 6 
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This fixes unique value of Jo = p e x,o u (L) for j = 1, 2 in the equation (jl.32p so that each <I>( J ) can 
be expressed as 

(1.41) z(A( Xn ) = l(J°)\[ P ' M ^ de =:H(j,V\x n » on [0,L]. 



2\pU)(x n )J Jk 

Since p"(p) > for p > 0, there exists unique p s > depending on Jo so that we have 

< for p < p s 



(1.42) 7J'(p) = l(pV(p)-Jo 2 ) 



= at p = p s 
> for p > p s 



It follow from (fP2j) that 



/ /J'(tp (1) + (1 -t)pW)dt > on [Q,L] 
Jo 



because p&> > p s on [Q,L] for j = 1,2 by 051) . Using (fTTTSl) . ([LT9D and ([OB with $(L) = =%o 
and Q ex = 0, one can easily show that $W — <I>( 2 ) satisfies 

($W _ $(2))" * * = on (o L), 

and the boundary conditions — $( 2 ) = at x n = and L. Then the maximum principle for 
elliptic equations implies = <I>( 2 ) on [0, L], and this yields p^ = p^ on [0, L]. So we get 
tpO) =<p&) on [0,L]. 

For any given (<5 en) o, ^o,o>£>ea;,o) £ ^Po> we have shown that there exists a unique one-dimensional 
subsonic solution (93,$). Thus, for the density p defined by (|1.14p from (</?,$), we get (pq,Eq) = 
(p($( 0), |^'(0)| 2 ), $'(0)) G qJi(Jo) where J is uniquely determined from ($ e n,o, ^0,0, Pexfi) through 
(jl.40p . Hence we can find the maximal set of parameter set ^3o satisfying the statement of Proposi- 
tion [L2] by collecting all parameters ($ en ,o> ^0,0 ,Pex,o) corresponding to (po, Eq) G U ^Pi(Jo)- □ 

Jo>0 

Remark 1.5. The proof of Proposition HOI shows that, for any given (<& en ,0i ^o,o,Pex,o) G *Po> the 
constant Jq(= pu) is uniquely determined. Therefore, one can use the first two equations in (jl.3ip 
to rewrite the last equation in (ll.3ip as 

(1.43) fa(& +p (p)j \ = p - bo . 

By (jl.34p . one can a/so find a unique constant p en ,o > satisfying the equation (jl.35p with po = 
Pen,o- Therefore, the boundary value problem in Proposition can be considered as a boundary 
value problem for the equation (|1.43j) with the boundary conditions 

P(0) = p enfl , p(L) = p exfi . 

And, this is equivalent to fix the pressure at both ends of the nozzle JV by (jl.3p . 

The rest of the paper is organized as follows. In Section we rewrite (|1.24p - (jl.27j) as a nonlinear 
boundary value problem for (ip, ^) := {ip — ifQ, <I> — &o), and set an iteration scheme to solve the 
nonlinear boundary value problem for (ip, ty). In Section O we achieve a priori H 1 estimate of weak 
solutions to a linear boundary value problem containing a second order elliptic system. By using 
this estimate, we prove the unique existence of weak solutions and establish a priori C a estimate 
of the weak solutions up to the boundary of M. Then a priori C 1,Q estimates and weighted C 2,a 
estimates of the weak solutions are obtained, and these estimates lead to prove the unique existence 
of C 2 solutions to the linear boundary value problem. It should be emphasized that the H 1 estimate 
of weak solutions for second order elliptic system is the key ingredient to get the main result of 
this paper. Section 0] is devoted to prove Theorem Q] by applying the a priori estimates given in 
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Section [3] and the contraction mapping principle. In Appendix [Aj the existence of one dimensional 
subsonic solutions to the Euler-Poisson system with variable background charge is proved. The last 
appendix gives an outline of proof for stability of subsonic flows under small perturbations of the 
boundaries of the nozzle N . 

2. Linearized boundary value problem and Iteration scheme 

For the rest of paper, let M, bo, L be as in TheoremCQ And we fix ($ en> Q, ^o,o>Pex,o) £ *Poi and 
let ((/9o,$o) be the corresponding background solution. Let b, (& en ,<& ex ,p ex ) satisfy the estimates 
(jl.28p for a G (0, a) with a to be determined later. 

2.1. Linearization of the equations in (jl.24p . For (z,q) = (z,qi,--- ,q n ) elxl™ and for p 
defined by (fTT4"l) . set 

(2.1) A(z,q) = {A x ,--- ,A n )(z,q) = p(z, |q| 2 )q, B(z, q) = p(z, |q| 2 ). 

Fix a constant Eo > 0. For any (z, q) satisfying z — ^|q| 2 > h(eo) for h defined by (II. lip . A(z, q) and 
B(z,q) in (I2.ip are well-defined, and they are differentiable with respect to z and q. Furthermore, 
one can directly check that 

B(z,q) 



d z Ai(z, q) = d z B(z, q)^, d z B(z, q) = 

(2.2) %A(2 . q) = B(2 . q) (^__||_ 

B{z,q)qj 



P'(B(z,q)) 



d q .B(z,c\) 



for i, j = 1, ■ ■ ■ ,n, where 5{j = 1 for i = j and 5y = otherwise. 

Suppose that ((p,&) G [C 2 (A/")] 2 satisfy the equations in (ll.24p in jV, and set 

(2.3) (^):=(^$)-(^o,$o). 

Since (</?q,$o) satisfy (jl.24p as well, by subtracting the equation div(p(3>0) |V(/2o| 2 )Vv?o) = from 
the first equation in (|1.24p . one has 

Lx(i>, tf) :=div I £ <9^($ , ityoW + *9 Z A($ , Dipo) 
( 2A ) \ j= i 

=divF(x,* ;J DV) in A/" 
for F = (Fi, • • • , F n ) given by 

(2.5) - Fi(x, z, q) = f z [d z M*o + ~z, D m + q)] ^L (0 ,o) + % [$a^i(*o, £>¥>o + q)] ~L ^ 

J 

where [G(X)] b x=a := G{b) — G{a) with <3?o and Dtpo evaluated at x G M. Here, D denotes 
(d xi , • • • , d Xn ), and each dj denotes d x . for j = 1, ■ ■ ■ , n. 

Next, subtracting the equation A$o = /K^Cb IV920I 2 ) — bo from the second equation in (ll.24p 
gives 

(2.6) L 2 (rp, *) := A* - 0,B($ O , D*^)* - 3 qJ B($ , £>¥>o) • £>V> = /(x, D^) + b - 6(x) 
for / given by 

(2.7) /(x, Zj q) = jT z [d z B($ + ~z, D m + q)] (S)L(o,o) + H [ d V B (®0, Dipo + q)] ^df. 

There is a constant <5i G (0,1) depending on the data such that each Fj(x, z,q) and /(x, z, q) are 
well defined for all (x, z, q) € for 

(2.8) P 35l := j\f x {(z,q) G R x R n : [z| + |q| < 3*i}. 
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Then, the following lemma can be easily checked. 

Lemma 2.1. Let Fi and f be defined by (|2.5p and (|2.7p . respectively. Then, 

(a) there exists a constant C depending only on the data such that 

F(x,0,0) = 0, /(x,0,0)=0, 

\D {z ^Fi{x,z^)\ < C{\z\ + [q|), |D ( , )q) /(x,2f,q)| < C{\z\ + |q|) 

for all (x, z, q) G TJ 2 s 1 and i = 1, • • • , n; 

(b) for any integer k > 2, there exists a constant Ck depending only on the data and k such 
that 

\D^Fi(x, z,q)| < C k , \D\ z ^f(x, z, q)| < C k 
for all (x, z, q) E T> 2 s 1 • Here, D*? z , denotes the k-th order derivatives with respect to z and 

q 

Lemma 2.2. Let A(z,q) and B(z,q) be as in (12. ip . 

(a) T/ie matrix [o>ij]ij = i := [<9 gj . Aj(<I>o, -tVo)]f?=i a strictly positive diagonal matrix in N, 
and there exits a constant A > satisfying 

1 

\' 

and such A depends only on the data. 

(b) Each aij is smooth in M . Also, there exists a constant C k > depending on the data and 
k to satisfy \\<Hj\\ckQ7) < for all i, j = 1, • • • , n and any nonnegative integer k. 

(c) Let h(p) be defined as in (jl.lip . For any constant eq > and any (z, q) £lx lR n satisfying 
z ~ IIqI 2 > h(e(j), the equality 

(2.10) 0*A(z,q) + q fl(*,q) = O 



(2-9) A/ n < [a^- =1 < T 4, 



(2.11) 



Proof. From (jl.38p and (|2.2p . it is clear that [ajj]" J=1 is a diagonal matrix with 

{B(® , D(p ) for i < n 

Furthermore, (|1.23p implies that there exists a constant A > satisfying (|2.9p . So (a) is proved, 
(b) is already stated in Definition 11.41 Finally, (c) can be easily checked from (12. 2p . □ 

Lemma [2.2f c) is crucial to get H 1 estimates of weak solutions to a linear boundary value problem 
associated with two elliptic operators L\ and L 2 defined by (|2.4p and (|2.6p . respectively. 



2.2. Boundary conditions for (vp, \E f ). Given <3? en , & ex ,p ex and ,^q, suppose that (ip, <£) satisfy 
the boundary conditions f|1.25|) (|l .2Tj) . Then (ip, ^) defined by (|2.3p satisfy 

ijj = on Tq, 



(2.12) 




- =%o) + ($en - $en,o) =: ^ en On T 
- ^0,o) + ®ex =■ ^ex On Y L 

d n v = o on r„, 



(2.13) p(B($o + y,V<po + V%l>))-p(B($o,V<po))=p ex -p exfi on V 

for B(z,q) defined by (|2.ip . To simplify notations, set 

p:= B{$ + y,V<p + Vi>), Pbg :=5($ ,V^o)- 
There exists a constant 5 2 > depending only on the data so that if 

(2-14) ll*llo°c*o + ^ 



ex 
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then 

(2.15) B($ + ¥, V(^ + i/>)) = h-\<S> + - i|V(v? + i>)\ 2 ) > > 

holds in M for the constant vq in ()1.23|) . Then the left-hand side of (|2.13[) is rewritten as (p — 
Pbg) Jo P'ftP + (1 — t)Pbg)dt, and p — pb g can be written as 
(P ~ Pbg) 

= * f B z ($ + **, V^o + tVip)dt + W • / £q($o + i*,V^ + *W>)cft. 
jo io 

It follows that the boundary condition (I2.13P can be rewritten as 

(2.16) fl q ($ , Vyo) ■ = Pex ~ Pex \ + g 2 (x, tt, W>) on T L 

for 

/■l 

(2.17) V^) = - y ^(^o + ^.V^ + tV^ + V^-tSq^o + ^V^ + iqJjjJ^ojttt. 

By (|2.12p . we can substitue \I/ = \I/ e:2: into (|2.16p to get a nonlinear boundary condition for ip: 

(2.18) S q ($ ,V^o)-VV = 5(x,VV,Pex,* e x) on T L 
for 

/ ,T, \ VeX ~ PeX,0 . A / yjy \ 

5(x,q,fe,*ex) = —7 — — +5f 2 (x,* ea; ,q) 

, * 5l(x,q,p ex ,* ex ) 
(2.iyj j 

5l (x,q,p ex ,* ex ) = / p / (t J B($o + ^ex,V<po + q) + (l-i) J B($o,V<y3o))c?t, 
Jo 

where &o,ipo, and are evaluated at x £ T^. 

Lemma 2.3. Let 5 2 be as in (I2.14|) . Under the assumption (jl.28j) . z/4<r < (52, £/ien i/te following 
properties hold: 

(a) £/iere exists a constant C depending only on the data so that 

|D q 5(x,q,p e:!; ,* ea; )| < C(\p ex (x) -p ex ,o\ + |*ex(x)| + |q|) 

for all (x, q) G x {q € K n : |q| < 25 2 }; 

(b) /or any integer k > 2, there exists a constant Ck depending only on the data and k so that 

\D*g(x,ci, Pex ,y ex )\ < C k 
for all (x, q) 6 Tl x {q £ R" : |q| < 25 2 }. 

Lemma 12.31 is a direct consequence of (12.15P and fj2 .191) so we skip the proof. 

For later use, we note that the boundary conditions d nw ip = on and f|2. 18p on T can be 
rewritten as conormal boundary conditions. On T w , the inward unit normal vector field n w on T w 
satisfies n w ■ (0, • • • , 0, 1) = 0. So, by (|2.1ip . d aw ip = is equivalent to 

(n-l) 

n 

(2.20) C^^jdjiP) -n w = on T w . 



It follows from (JZ2D that 5 q ($ , Vy> ) = (0, ■ ■ ■ ,0, $ J o ° Vyo)) ) for J = B(* 0) V^ )d„^ where 

(n-l) 

^o^o and <9 n (/Jo are evaluated on T^. Combining this with (|2.18p gives 

( 2 - 21 ) O^o-ijOj^) ■ n L = g(x, V^,p ex ,^ ex ) on T L 

., Jo 
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for the inward unit normal = (0, • • • ,0, —1) on Tl 




2.3. Iteration scheme. Suppose that a < 1 in Theorem [TJ In £12.11 and §2.21 we have seen that 
(tp, $) G [C^AO n C 2 {M)} 2 solve (0311 - 071) if and only if (ip, *) = (<p - ip , $ - $ ) satisfy the 
equations f|2.4j) and (|2.6p in A/" and the boundary conditions (|2.12p and f|2. 18j) provided that 

(2-22) ll*llc°(AD + ll^llci(AT) ^ 2min{<Ji,*2} 

for 5\ and (52 in f|2.8j) and (|2.14p . respectively. 

Given q G (0, 1), we define an iteration set Km as follows: 

(2.23) /C M := K$ x 4? 
where /C^ and /C^ 2 ) are given by 

Kff := € C^ 1 _ <> , rouri) M : * = on r„, ||^ ||^ 1 JJ^« Uri > < f <r} 

*S? := {* € C£_ arouri) (A0 : ||*||£-«™*'> < £ ,> 

for constants M > 1 and tr G (0, 1) to be determined later. If Ma < 2 min{o~i, o~2}, then F(x, ^, Df/"), 
/(x, \E', D?/>) are well defined in M for all (ip, *) G Km, and g(x, V?/ 1 iVexi ^ ex) is well defined on Tl 
for F, /, g defined by (|2.5p . (|2.7p . (|2.19p . respectively. Consider the following linear boundary value 
problem for (if), ^f) 

(2.24) fL^,*) = divF(x,* L ^) 

with boundary conditions (12. 12[) (replacing (ip^) 011 the left-hand side by (ip,*&)) and 

(2.25) £ q ($o,V^o)-V^ = 5 (x,V^,fe,*ez) on T L , 

where L\ and L2 are defined in (I2.4p and ()2.6p . respectively. If the linear elliptic system (I2.24D with 
boundary conditions (|2.12p and (|2.25p has a unique solution G [C^AT) n C 2 (A/")] 2 , then an 

iteration mapping J can be defined by 

(2.26) 3: K M ) 

It is easy to see that (tp* , G /Cm is a fixed point of J if and only if (ip* , $*) = (930, ^0) + (V'*, ^*) 
is a solution to the nonlinear boundary value problem of (|1.24p - (|1.27p . So it suffices to prove 
unique existence of fixed point of 3 in /Cm for the proof of Theorem [TJ 

Proposition 2.4. Let J\f be as in (jl.6p . Fix 60 > and L > 0, and let the parameter set be 
as in Proposition Given ($ e n,o> ^o,o,Pex,o) G tyo, let (<£>o>^o) ^ e ^ e corresponding background 
solution. Assume that b, (<& en ,& ex ,p ex ) and satisfy (|1.28p and (jl.29j) . T/ien, /or any given 
a G (0,1), i/iere exisi constants M > 1 and ai G (0,1) depending on the data and a such that 
wherever a G (0, <n] in (|1.28p . t/ie iteration mapping 3 defined by ()2.26j) /ias a unique (tp*,^f*) in 
Km satisfying 

3(ip*,^*) = (tp*,V*), 

(2.27) p>($ + |V^o + W*| 2 )) - |V^ + W*| 2 > vx 
for a constant v\ > 0. 

Theorem [TJ easily follows from Proposition 12.41 by choosing = o~\. So the rest of the paper is 
devoted to prove Proposition 12.41 



SUBSONIC FLOW FOR MULTIDIMENSIONAL EULER-POISSON SYSTEM 



13 



3. Gradient estimates for elliptic system and mixed boundary value problem 
The following a priori estimate is essential to prove Proposition 12.41 

Proposition 3.1. Fix a G (0,1) and functions F G C}f a ToUTL) (M',^ n ), f G C a (N), W en G 
C 2 ' a (ro), W ex G C 2 ' a (r^) and a function g G C a (r^) satisfying the following properties: 

(3.1) F • 6j = on r /or j = 1, • • • , n - 1; 

(3.2) l|F||iXjJ oUri) + H/lUx + IblU.FL + IIW^H Wo + \\We X h,a,r L < 1 
where ej denotes the constant vector (0, • • • , 0,^1^,, 0, • • • ,0) for each j = 1, • • ■ , n — 1; 

(3.3) d nw w en = o on to n 9 nro w ea; = o on rT n iv 

Then the linear boundary value problem 

(3.4) |W) = divF 

(3.5) » = on T , 5 nil ,i' = on T w , B q ($o, V^) • Vt) = 3 on T L 

(3.6) W=lZ en "I"' 9 Dw W = 0onT w 

\W ex on Tl 

has a unique solution (v,W) G [C^!°[_ a r ur L )(-^)] 2 satisfying the estimate 

(3.7) \H { C^ T0VJVl) + IWfc 1 ^ ^ < C*, 
where the constant depends only on the data and a. 

Corollary 3.2. // the condition (|3.2p is dropped in Proposition I3.il i/ien t/ie linear boundary 
value problem of ()3.4|) f)3.6|) /ias a unique solution (v,W) G [C^"_ Q r ur L )(-^)] 2 ™^ satisfying 
the estimate 



Proof. Set m := ||F||J ^° UFl) + \\f\\aM+ II«IUtl + ll^nllwo + HWexlIWi- Then ^ follows 



I, |,(-i-a,r ur L ) n wr ||(-i-a ) r ur i ) 

^ 38 ^ W V h, a ,M + W W h,a,jV 

< CH\\n ( d 0UrL) + WfWaM + \\9\\a,T L + ||W Bn || 2 ,a,ro + \\W ex \\ 2 ,a,r L ) 

for the constant C* same as in Proposition \3.1[ 
l(-a,r, 

by applying Proposition 13 . 1 1 to ^). □ 

In order to prove Proposition 13. H first of all, we achieve uniform if 1 estimate of (t> , W) under the 
assumption of (|3.2p . then we combine the uniform H 1 estimate with 1? integral growth estimate of 
(Dv, DW) to get uniform C a estimate of (y, W) in N . H 1 and C a estimate of (v, W) are obtained 
through estimates of weak solutions to the linear elliptic system ()3.4p . For uniform estimate of 
(v, W) in C 1,Q -norms or higher, we apply a priori C a estimate of weak solutions and use individual 
elliptic equations of the system (|3.4p in the correct order. 

3.1. H 1 estimates. Define 

Ui := {£ G H\Af) : ^ = on r }, U 2 ■= {v € B x (tf) : 77 = on T U T L }. 
Suppose that (v, W) G [C l (N) n C 2 {M)} 2 solve (|5^j> - (|33|l then, for any (£, 77) G Hi x we have 



(3.9) 



/ Li(i7,W)£ + -&2(«,W>7<fc = / (divF)^ + fv dx.. 
JN JjV 
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For the fixed functions W en and W ex in (|3.6p , we define a function W bd in M by 
(3.10) W bd (x) = (1 - ^)W en (x') + ^W^x') 

where x' € A for x = (x', x n ) € TV. Then, for any (£, 7?) G %i x H 2 , (v, W) := (v, W) - (0, W M ) G 
Hi x 7{ 2 satisfy 

(3.H) C[(v, W), (£, r,)} = <(F, /, g, W bd ), (£, ?7)) 

where 

£[(«, W), (£, r?)] := A[(t5, W), £] + C 2 [(v, W),rj], 
(3.12) := J aud&dtt + WdtA^Difti-DZdx, 

£z[(v,W),rj\ := / W-,Dr ? + (#9 2 S(^ ,^o) + 9 q S(* ,^ )-^)r ? dx, 



and 



((F, /, g, W bd ), & ??)) := ((F) /, s), (£, ^)}l + <W M , (£, r?)} ; 



(=:/(!)) (=:/( 2 )) 

(B($ ,V^))^ 



dS 



(3-13) /W(C,r ? )= [ F-Dt-fr,)**.- [ (F-a^+ / W (g(*o, 
I (2) (Z,V) = - I W bd d z A(* ,D m ) • Dt + (W bd d z B(y ,Dp ) - &W bd )r, dx 

•AAA 

for unit outward normal vector n OU £ on 97V \ To- This can be directly checked by applying the 
integration by parts to (|3.9p and using (|2.2ip . In (|3.12p and hereafter, the Einstein summation 
convention is used. We note that the compatibility condition (|3.3p implies d nw W bd = on T w . 



Definition 3.3. We call (v,W) a weak solution to (ROD -rtTTTH) if (v,W) = (v,W) - (0,W bd ) G 
Hi x H 2 safe/?/ dHU]) /or any (£,77) e^x%. 

Lemma 3.4. Under the same assumption as Proposition \3.1l the linear boundary value problem 
(pTi]) - (j3"l)} /ias a unioue wea& solution (v,W) G [id^A/")] 2 , and (5,1^) := (u,W) - (0,W M ) satisfy 



(3-14) IHIffi(A0 + II^IUi(AA) < Ch 

for a constant Ch depending only on the data. 

Proof. Note that H := 7di x H 2 is a Hilbert space equipped with an inner product (•, defined 
by 

((Ci,m),(^m))H= I (D^-D^ + Drjt-D^dx. 
JN 

By Poincare inequality, there exists a constant (3q > satisfying 
(3-15) A)(lieil2rxcA0 + Nl 

for any (£, n) € H, and such a constant [3q can be chosen depending only on n and A/". 

Thanks to Lemma l2T27c) . we have r]d z A{^ ,Dtp ) ■ D£ + r]dqB(^f , Dip ) ■ D£ = for any 
(£, 77) G H. Then the definition of £[•, •] in (I3T2H yields 

(3.16) £[(£, 7/), (£, r?)] = / a^c^ + l^l 2 + 0*B(* O , £<A))?7 2 dx. 

JN 

Combining (12.2p . Lemma 12.2( a) and (11.23P together, we also get 

(3.17) £[(£,77), (£,77)] >A / |I^| 2 + |^| 2 dx 
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for Ao = min(A, 1) where A is the one appeared in (|2.9p . On the right-hand side of (|3.1ip . applying 
Holder inequality, trace inequality and Poincare inequality to ((F, /, g, Wm), one can find a 

constant C n h > depending only on n and AT to satisfy the estimates 

1/2 

(3.18) 



J |ZA£| 2 + \D V \ 2 dx^j 



This can be directly checked from (|3.13|) . We note that C n h is independent of F,f,g and Wm 
because of (13. 2p and (I3,10p . Then, the Lax-Milgram theorem implies that there exists unique 
(v,W) G V. satisfying (fBTTTT) for all (£,r/) G U. Finally, substituting (f,r?) = (v,W) into (133x1) . and 
applying (|3.17p and (|3. 18[) give 

1/2 

(3.19) A 



\Dv\ 2 + \DW\ 2 dx < C nh 



AT 



\Dv\ 2 + \DW\ 2 dx 



N 



(|3.14p follows from applying Holder inequality and Poincare inequality to (|3.19p . 
3.2. C a estimate. 



□ 



Lemma 3.5. Under the same assumption as Proposition \3. 11 let (v, W) be as in Lemma \3.4\ Then, 
for any a G (0,1), there exists a constant Ca depending only on the data and a so that (v,W) 
satisfy 

(3.20) \\v\\a,M+\\W\\ a ,M <Ca- 

Proof. Since Af is a cylindrical domain with the cross-section A, there is a constant Rq > depend- 
ing only on n and A such that for any x* S M and r € (0, Ro], the inequality 

i_ < vol(ff r (x,)nAQ < Cq 
Co ~ vol(B r (x*)) 

holds for a uniform constant Co > 0. So if there are constants R\ < Rq and kq > satisfying 



(3.21) 



Dv\ 2 + \DW\ Z dx < Kffr 



2n-2+2a 



'B r (x»)rW 

for all x* G M and r G (0, R\], then it follows from |13|. Theorem 3.1] that 

\v\a,M + \W\ a ,M <C \K + \\v\\l2(AT) + HW'llia^ 
for a constant C depending only on n,a, A and L. Combining this with Lemma 13.41 gives 

\v\u,N + \W\ a ,N < C(k + C H ), 

and this estimate proves (|3.20p . To obtain (|3.2ip . we need to consider three cases: (i) B r (x) C M, 
(ii) x G dN \ C for C := (To" U Tl) n F^, and B r (x) n C = 0, (iii) x G C. More general cases can 
be treated via these three cases. Also, cases (i) and (ii) are easier to handle than case (iii). So we 
treat only case (iii) here. 

Fix xo G C. Without loss of generality, we assume xo G C n Tl. For a constant r with < r < 
jq min(L, diam A, 1), set T> r (xo) := M Pi B t (xq). Let w\ and ^2 be weak solutions of the following 
problems: 

di(J2]=i a ij( x o)djWi) = in P r (x ), 
(YTj=i aij(x )djWi) ■ n out = on <9r> r (x ) n (T L U F w ), 
w\ = v on dV r (xo) nTV, 



and 



Aw2 = in P r (xo), 

dn out w 2 = o on ap r (x ) n (r L u r™), 

w 2 = W on «9P r (x ) HA/", 
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for unit outward normal vector on dV r fl (r^ U T w ). In other words, (101,102) satisfy 

(3.22) / a,ij(xo)djWidiZi + D1V2 ■ Dz 2 dx = 

•/XV(xo) 

for any z\,z% £ H l {T> r {x$)) n {z = on dT> T (xo) fl M}. We extend z\ and 22 by setting Z\ = Z2 = 
in jV\2?,.(a;o), and substitute = (#1,22) into (13. lip . Then subtracting (I3.22p from (13. lip with 

(£,77) = (21,22) gives 



/ aij(x )di(v — w\)djZ\ + D(W — w 2 ) ■ Dz 2 dx 
(3.23) •m-(xo) 

= <(F,/, & ,TV M ) J (zi,^)>- / Tdx, 

JV r (x ) 



where T is defined by T = 71 + T 2 + T3 with 

71 = (ay(xo) - aij(x))djvdiZx, T2 = Wd z A(^ , Dip ) ■ Dz 1} 

r 3 = (iya ZJ B(^ , Apo) + q £(* , £Vo) • ihj)^. 

Set (Vi, V2) := (v,W) - (wi,w 2 ). We substitute (21,22) = (Ul, V 2 ) into (ET231 and use ([23]) to get 

(3.24) / aij (xoJflMfyVi + DV 2 ■ DV 2 dx>\ I \DV l \ 2 + \DV 2 \ 2 dx 

JT> r (x ) JV r (x ) 

for Ao same as in (|3.17p . It remains to estimate the right-hand side of (|3.23p . Using Lemma l2.2f b) 
and Holder inequality gives 

/ \ 1/2 / v 1/2 

(3.25) / \Ti\ dx< |K|| c insn K / \Dv\ 2 dx) If |Wi| 2 dx) . 



v r (xo) \ Jv r (x ) J yjx> r (x ) 

Similarly, the following inequalities are obtained by Holder inequality and Poincare inequality: 

/ \T 2 \ dx<in (r 2 f \DW\ 2 dx) if |Wi| 2 dx] 
(3.26) \ 1/2 1/2 

\T 3 \ dx < 4/xi I r 2 f \Dv\ 2 + \DW\ 2 dx J [/ |/JU 2 | 2 dxj 



^(lo) \ JT> r (x ) J \JT>r(x ) 

for /ii = ||S*A(* ,-D^)IU«(A0 + l|fl*-B(*o,^)lk»(AO + K^o, £Vo)||l~(aO- From ([535]), 
the estimate |((F, /,<?, W M , (Fi, F 2 )))| < |/W(Vi, U 2 )| + |/ {2) (^i, ^2) I follows, and \I^{V U V 2 )\ < 
\J\ \ + holds where 

Ji= / (F • DV\ — fV 2 ) dx, 
JT> r {xo) 

JdVr(x )ndM Jdv r (xo)nr L J o 

By using Holder inequality, (|3.2p and Lemma l2.2l fb). we get 

(3.27) \J2\<Cr^(f \V!\ 2 dx\ , 

\ia© r (x ) y 

where the constant C in (|3.27p depends only on the data. Unless otherwise specified, we presume 
that any constant C appearing in various estimates depends only on the data for the rest of the 
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paper even though C may be different in each estimate. Since V\ = on dT> r (xo) DJ\f, one can 
apply the trace inequality and Poincare inequality with scaling to obtain 

/ \ 1/2 / \ 1/2 

/ iV^dS) <C[r[ |Wi| 2 <fc) , 

\Jav r { XQ ) J \ Jv r (-K ) J 

and this combined with (I3.27P gives 

\ 1/2 



2 dx 



(3.28) \J 2 \ < Cr n l 2 [ [ |Wi| 

Vyxv(xo) 

Also, it is easy to show 

/ \ 1/2 

(3.29) |Ji| <Cr n l 2 (f |Wi| 2 + |W 2 | 2 dxj . 

\J-D r (xo) J 

By ([330]) and ([321), W 6d satisfies 

(3.30) IIW'Mlb.a.JV < C, 
and this yields 

/ \ 1/2 

(3.31) \I^(V U V 2 )\ < Cr n ' 2 [ [ [Wil 2 + \DV 2 \ 2 dx) 

\JVr(xo) J 

It follows from (^2i]l - (13361) . (13381) . (13391) and (133TD that 



\DVi\ J + \DV 2 \ Z dx<C\r z / \Dv\ z + |Z)VT| 2 + r z 



then |13l Corollary 3.11] implies that 



\Dv? + \DW\' Z dx < C I ( — + r 2 ) I \Dv\ 2 + \DW\ 2 dx + r n 



'Cp(xo) \ V / JVr(x ) ) 

for any < p < r(< min(L, diamA, 1)). Then, by |13} Lemma 3.4], for any a 6 (0,1), one 
can find a constant tq with < rg < ^min(L,diam A, 1) depending on the data and a so that 
wherever < r < ro, (v, W) satisfy 

(3.32) f \Dv\ 2 + \DW\ 2 dx < Cr n ~ 2+2a for any x £CnT[. 

JV r (xo) 

In fact, one can similarly show that (I3.32p holds for all xq E N . □ 

3.3. C > a estimates. Under the assumption of Proposition 13. 1\ let (v, W) be a weak solution of 
(|3.4p - (|3.6p in the sense of Definition 13.31 Then Lemma 13.51 combined with (j3.30p provides C a 
estimate of (v , W) = (v, W) + (0, Wm). In order to get uniform C 1,a estimates of (v, W), we regard 
the elliptic system (13.4p as two separate elliptic equations. This is possible because the principal 
parts of (j3.4p are not coupled. The Holder gradient estimate for first equation in (j3.4p gives uniform 
C ,a estimate of v. Then C ' a estimate of W is obtained by applying Lemma [3.51 and uniform C > a 
estimate of v. 

Lemma 3.6. Under the same assumption as Proposition \3.1[ let (v, W) be as in Lemma \3.4\ Then, 
for any a 6 (0,1), there exists a constant Cb depending only on the data and a so that (v,W) 
satisfy 

(3.33) |b||i,aJV+||W||i, a jv<C7B. 
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Proof. By (pOOl) . it suffices to show that {v, W) = (v, W) - (0, W bd ) satisfy 
(3.34) \\v\\i,a^ + \\W\\i, a ^ < C 

for a constant C depending only on the data and a. 

Step 1. By choosing r] = in (13,111) . v can be regarded as a weak solution to 

'$(£i=i HA*) = div ( F " Wd z A($ , Dipo)) =: divF* in TV, 



(3.35) 

By Lemma 
(3.36) 



v = 

(Z)i=i a ijdjv) ■ n w = 



TO 



F* satisfies the estimate 



annp'{B(3> ,Vip ))g , 

Jo ~- 



on r , 

on T w , 
on T L . 



\F*\Lm<C. 



Since f) G C a (r^), one can find a function Vh G C ,a (AAn {x n > -j}) satisfying 



(3.37) 



(y^QjjQj^) • n L = f) on F L and 



l,a^Vn{a;„>f } 



< c. 



One can refer to Step 2 in the proof of Lemma 3.1 of [3] for more details on the existence of such 
vi,. It is easy to see that Vj, := v — v§ is a weak solution to 



(3.38) 



f fle(£"=i Otf W = div ( F * " E"=i Oij^) =: divH in TVn {x n > L/4}, 

(J2j=i a ijdjVt>) ■ n w = -QZ"=i a ij d 3 v t)) ' n «> =: si ° n r «> n {x n > L/4}, 

Vf, = v - Vf, =: g 2 on AT n {x n = L/4}, 

> (E"=i a y 9 j H) • = on T L . 

For | € K + , set Q+ := A x (Z,2L - /). For (x',x„) G Qt, define 



(3.39) 

for i,j = 1, 
in Qt by 



ajj(x',x n ) = chij(x',2L - x n ), g k {x',x n ) = g k {x',2L - x n ) 
, n, k = 1, 2 and for all x' = (xi, ■ ■ ■ , x n —\) G A. Also, we define H = {Hi, ■ ■ ■ , H n ) 



(3.40) 

and 

(3-41) 



for x n < L 



TT i i \ I Hj{x , x n ) .... .. ,, _ 

HAx,x n ) = < ~ , , „ for 7 7^ n, 

JV 7 \fT i (x , 1 2Z/-x„) forx n >L 



H n (x , x n ) 



H n {x', x n ) - H n (x', L) for x n < L, 

-{H n {x f , 2L - x n ) - # n (x', L)) for x n > L. 



According to ()3.39p and (|3.40p . aij,gk and each Hj(j ^ n) are extended from A/"n{x n > j} to Q 



by even reflection about x n = L, and H n is defined as the odd extension of H n {x', x n ) — H n {x', L) 

about x n = L. Hence ay G C ' 1 ^ ), H G C a (Ql) and 51 G C a (<9A x {|x n - L| < 3L/4}). For 

4 4 

such extended functions 0^,51,52 and H, the boundary value problem 

'^(E"=i%^^) = divH, in Q+ 
< 3. 12 i <J (£» =1 ay^.f/) • n w = ^ on «9A x {|x n - L| < 3L/4}, 

U = c/2 on A x {|x n - L| = 3L/4} 
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has a unique weak solution U € H l (Q\). Furthermore, the assumption (j3.2j) together with the 
proof of [31 Lemma 3.1] shows that 

(3-43) W U k a> Ql <C(\\U\\ l2{q+l) + 1). 



T 



It also follows from ()3.39p and (|3.40p that if U is a weak solution of (|3.42p . then so is U(x',x n ) = 
U(x',2L — x n ). Then the uniqueness of weak solution yields U(x',x n ) = U(x',2L — x n ) thus 
d n U(x',L) = which implies (Y^=i a ij^j^) ' n L = on r^,. From this, Vf, = U is obtained in 
Mn{\x n -L\< 3L/4}. Therefore, by Lemma[3ja j33SD, (f3737|^ and (HEMP , the estimate 

holds for a constant C depending only on the data and a. 

Step 2. Back to (|53S)» . by flES), Definition O and (J3TD, F* satisfies 

(3.45) F* • % = on r for j = 1, ■ ■ ■ , n - 1. 

For / G M + , set := A x (—1,1), and we extend and F* into Q 3L u as follows: for x n < 0, 
define 

(3.46) aij(x,x n ) = aij(x', -x n ) for x'eA, i,j = 1, ■■■ ,n. 
Define F = (Fx, ■ ■ ■ ,F n ) in Q~ L/4 by 

£,// \ jF*{x',x n ) forx n >0 
-F 1 (x,x n ) = < J , , x tor i 7^ n, 

' l ' ^ \-F;(x',-x n ) for x n <0 J ^ 



(3.47) 



F n (x , x r 



F*(x',x n ) for x n > 
F*(x', -x n ) for x n < 



It follows from (|3.2p and (|3.45p that each Fj is in C a (Q 3L ^ 4 ) and satisfies \\Fj \\a,Q 3L/4 < C. Hence, 
the boundary value problem 



(3.48) 



diQ2]=i o-ijdjU) = div F in Q 3L/4 

(E?=i OiiW • = on9Ax^ ^ 



4 ' 4 



U(x',x n ) = (sgn a;„)£(a/, |x n |) on A x {\x n \ = ^} 
has a unique weak solution [/ E H 1 ^^^). Moreover, C/ satisfies the estimate 

(3-49) rilw^ < C. 

In fact, (|3.43p and (|3.49p can be verified by simplifying the estimates of [17J to the case of linear 
elliptic equation (see [T7| Theorems 1.1 and 1.2]). By using (|3.46p . f|3.47|) and the uniqueness of 
weak solution to ()3.48p . one can directly check that U(x',—x n ) = —U(x',x n ) in Q 3L u, and this 

yields U(x',0) = on T . Therefore, we get v = U in J\fn {x n < 3L/4}. It follows from (13301 and 
(|3.49p that ||u||i i a i A/'n{a; n <2L/3} < C. Combining this with (|3.44j) yields 

(3.50) \\v\\i,*M < C 

for a constant C depending only on the data. 

Step 3. Substituting £ = into (|3.1ip . W can be regarded as a weak solution to 



(3.51) 



= f + {W+W u )d z B{$ Q ,D<p Q )-mr u + d^B{$ ,Dy )-Dv in M, 

w = o, on r ur L , 

DW n w = 0, on T w . 
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By using odd extensions of f* and V about Tq and Tl, one can argue similarly as Step 2 to show- 
that 

(3.52) \\W\\ 1>aM < C 

for a constant C depending only on the data and a. □ 

3.4. Unique solvability of mixed boundary value problem. In £ j3.1| - £|3.3[ it is shown that the 
boundary value problem (j3.4|) ()3.6[) has a unique weak solution (v, W) and that the weak solution 
is C 1,a up to the boundary of Af. In order to prove Proposition 13.11 it remains to show that v and 
w are C 2 in M thus the weak solution is actually a classical solution of (|3.4j) — ()3.6|) . 



Proof of Proposition \3.1\ By Lemmas I3.4H3.61 the linear boundary value problem ()3.4j) — f)3.6|) has 
a unique weak solution (v, W) in the sense of Definition 13.31 and (v, W) satisfy the C ,a estimate 
(|3.33p . Then v can be regarded as a weak solution of 



di(^2 aijdjv) = div(F - Wd z A(<S> , Dip )) =: divF in M 
(3-53) 

C/ y aijdjv) ■ n w = on T w 

3=1 

for F G C 1,ct (A/' U T w ). Since every aij is smooth in Af, v is in C 2 ' a (Af) . In fact, v is C 2,a up 
to r„, away from Tq U Tl- We give a heuristic argument to verify this. For any given point 
xo G \ (ro U T^), can be locally flattened near xo, and (I3.53j) can be written as a conormal 
boundary value problem as follows: 



n 



(3-54) 



'V>><,'V ! = divfl in B 2R (°) ■■= B 2RoW n {y n > 0}, 

n 

C^2 ^ijdy^) • n w = on dB+ Ro (0) n {y n = 0} 



j-- 



for sufficiently small Ro > and inward unit normal vector h w on {y n = 0}. We note that h w 
is a constant vector. Here, we assume that xo is mapped to in the flattened domain. Since the 
boundary of the cross-section A of M is assumed to be smooth, one can find a smooth flattening 
map near xo so that dij G C 3 (B Ro (0)) and F G C 1,a (B Ro (0)). For each k = 1, • • • , n — 1, we 
differentiate (|3.54p with respect to y^ then w^ k \y) := v Vk {y) becomes a weak solution to 

n n 

fytCE ~ a ^vM k) ) = div^F - £ d yk ~a l3 d y] v) =: div.G^y) in 5+ (0) 
3=1 i=i 

n n 

E ^A^) • n w = dy^dy.v) -n w =: g^(y) on dB+ Q (0) n {y n = 0} 

3=1 3=1 



for G( x °) G C q (£+ o (0)) and <?( x °) G C^B^O)). Thus u;(*) G C 1 ' a (S+ )/2 (0)) for k = 1, ■ ■ ■ ,n-l, 
and this implies that is C a up to dB^ (0) unless i = j = n. One can also check that v ynVn is 
C a up to dBft o (0) by solving the first equation in (|3.54p for v ynVn because a nn is strictly positive 

in B2 Rq (0). This proves that v is C 2,a up to T w away from I^o U Tl. Rigorous proof can be given 
by using difference quotient (cf.JTTJ Section 8.3]). Once v G C 2,a is shown, we can also show that 
W is C 2 ' a up to away from To UTl using the second equation of (|3.4p and C 2,a regularity of v. 
It remains to estimate C 2 ^_ a p oUr ^-norms of v and W by a scaling argument. For a fixed point 
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yo G AT \ (ro UTl), let 2d = dist(yo, Tq UT/,), and set a scaled function 

^ (yo) («) == J+^[<y° + ds ) " "fro) - dDv(y ) ■ s) 

for s£{s£ -Bi(O) : yo + ds G A/"} =: fii(yo) where -Bi(O) is the unit ball in W 1 with the center at 
the origin. By (|3.33p . we have 

(3-55) lk (2/o) llc70( Bl (Q)) < \\v\\ 1:aM <C B . 

Substituting into (|3T53|> gives 

'^(Ey=i%-(yo + ^)5 s , v W) = div s F(yo) in fii(o), 



(3.56) 
for 

(3.57) 



(E;=i (yo + <k)d s .u) • n„,(y + ds) = fl (*>), on 0Bi(O) n {y + ds G 

[ (ifaB 1 (o)n{yo + d S Gr w }^0) 

p( yo )/ n ._ F(y + ds) - F(y ) - g " =1 (aij(yo + <*s) ~ aij(j ))djv(j ) 



d° 

i(yo) , (E"=i °ii(yo + ds)<9,-v(y )) • n^yo + (is) - (E"=i a ij{yo)djv{y )) ■ n w (y ) 



d a 

By using Lemma l2T2T b). (|3.2j) and (|3.33[) . one can find a constant Cp depending only on the data 
and a so that, for any y G Jf \ (T U T L ), F^ and g M satisfy 

(3-58) ||F^)|| liQ)Bl(0 ) + ||g (yo) ||i, Q , Bl( o) < C F . 

Then, it follows from [11, Sections 6.3 and 6.7], (1335]) and (1338]) that is in C 2 > a (B 1/2 (0)) and 
satisfies 

lb (yo) l| 2 , a , Bl/2 (o) <c(c B + c F ) =-.a 

for the constant Cb in ()3.33p . We emphasize that the constant C* above is independent of yo G 
JV \ (ro U r^). Therefore, re-scaling back gives 

II, ,||(-i-a,rour £ ) ^ r 

for a constant C** depending only on the data and a. 
Similarly, we consider W as a solution of the equation 

AW = f + d,B($ ,D<po)W + dqB($ ,Dtpo)>Dv=:U in M. 

Since f* G C a {M) with satisfying ||f*|| a) AT < C by (13.33j) . we can repeat the scaling argument above 
to conclude that 

wt,U' To[JVL) < c *- 

□ 



4. Proof of Proposition 12.41 

Fix a G (0,1). If Ma < 2min{<5i,5 2 } for 61,62 in (i!T8j) and (I2.14D respectively, then for each 
(■0, \P) G /Cm, functions F(x, ^, Dip), f(x,^,Dtp) are well defined in M, and c/(x, Dip,p ex , ^ e x) is 
well defined on T L where F(x, 2, q), /(x, z, q) and <?(x,q, 

Pexi^ ex) defined by (|2.5p . (|2.7p and 

(I2.19p . respectively. 

Lemma 4.1. There exists a positive constant ^3(< 1) depending only on the data so that if Ma < £3 
in the definition of K,m given in (I2.23p . i/ien i/ie following properties hold: for any (ip,^) G ICm, 

(4.1) ||F(x,¥, J^)||£# oUr,l) < C(Maf 

(4.2) ||/(x,^,^)|| QX <C(M 2 (7 + l)a 

(4.3) || 5 (x, Di>,p ex ,^ ex )\\ a ^ L <C(M 2 a + l)a 
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for a constant C depending only on the data. Furthermore, we have 

(4.4) F(x,#,£ty!) ■% = on T G for j = 1, • • • , n - 1. 

Proof. For /? defined by (|1.14|) . we can choose a constant J > depending only on the data so that 
if Ma < 5, then there holds 

0< -infp($ ,D<po) < p(<S>o + y,Dp + D4>) < 2supp($ ,-CVo) m M 
2 M j\f 

for all (tp,^) G /Cm- Then, by (|1.4p . the inequality 

Mo < p'(p($o + *, £>¥>o + £>V0) < — m M 

Mo 

holds for a constant /xq > depending only on the data. Thus (|4.ip - (|4.3p easily follow from (jl.28p . 
(1231) . (H2D, (I2TT6D . pTTjl . Lemma O and LemmallOJ 

In order to prove (|4.4p , we need to take a closer look at the definition of F in (|2.5p . By (|2.2p , we 
have c^A(<3?o + z, Dfo + q) || Dip® + q for all (x, z, q) € where the set P35 1 is defined by (|2.8p . 
Choosing 5 3 = min{<5, <5i, £2} so that for any (x, z, q) G AAx {(z, q) G MxR n : |z| + |q| < #3}(=: X^), 
we have 

I' 1 

j z[d z Ai($ + z, Dtp + q)](2^) = ( 0)O ) di = °i( x ' z > q)- ^ + 02 (x, z, q)q, 

/ S~] qj[d qj Ai($ , Dp + q)]~'L di = a 3)i (x, 2, q)aj for all i ^ n. 
^ f=l 



(4.5) 



for 01,02,03,1 G C(Vg 3 ,'K). For any (ip,^) G /Cm, since ^ = on To by the definition of /Cm in 
(j2T23|) . there holds 0^ = on T for all i^n. So if Ma < 5 3 , then, by fl23J and g3J, F(x, D^) 
satisfies glj. □ 

Now we are in position to prove Proposition 12.41 

Proof of Proposition \K4\ If M 2 a < 1, then, it follows from (ll.28p . Proposition 13. II and Lemma 14. II 
that the elliptic system (|2.24p with boundary conditions of (|2.12p and (|2.25p has a unique solution 
G [C^.^ur^N)] 2 , which satisfies 

(4.6) Mt^'^ + ||*||^ roUr£) < 2&a 

for C* depending only on the data and a. 

Next, we estimate - J(^2,f 2) in C^i^rour^)^) for e = M)- For 

each fc = 1, 2, set 

(4,*ifc) := 3($k,*k) in AT 

F W := F(x, ¥ fc , £ty fc ), /« := /(x, tf fc , £ty fc ) in TV 
o (fc) := g(x,Dtp k ,p ex ,^ ex ) on r ex . 

Then (u*,iy*) := $1,4! 1) - (^2,^2) satisfy 

JiiCu*,^*) = div(F« -F( 2 )) 

\l 2 (,*,^) = /( 1 )-/( 2 ) 

w* = on r , B n ($ ,D<po) -Dv* = - 5 (2) on T L , 
dn„y = 3n„W* = on r WJ 

w* = o on r ur L . 
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By Lemma [2jja) and Lemma [2T3l f a) . there exists a constant C depending on the data and a such 
that 

l|F« -F( 2 )||^° Uri) + -/ (2 W+ \\ 9 U -9^\\ a ,r L 

47 < C(l + M)a{U, - U { C7 ,T ^ VL) + H*i - ^2||& Q ' r ° Uri) )- 

By Corollary 13. 2l and (|4.7j) . one can find a constant depending on the data and a so that (v* , W*) 
satisfy 

h *ii(-i-o,r ur Zv ) n w *n(-i-a,rour L ) 

< C\l + M M ||^ - ^2||^' r0Uri) + 11*1 - *2||fc 1 7 ,r ° UrL) )- 

To complete the proof, we choose M and a as follows: 

(i) Choose M to be 

(4.9) M = max{l,4C b } 

for as in (|4.6p . Then 3 maps /Cm into itself. 

(ii) Next we choose u\ as 

(4.10) u ' = m^ n{hM ^W2Mc^ 

for 81,82, 83 as in (|2.8p , (|2.14p , Lemma 14. 1^ respectively. Then (|4,8p implies that 3 is 
contracting in [C/ 1 _ a r urx,) 

(A/")] 2 - 

For M as in (14.9P and for any o~ G (0, o~i], /Cm is a nonempty Banach space equipped with the norm 
°f II ' II2 Q 1 A^' r ° Uri ^ Then the contraction mapping principle implies that for any a G (0, o"i], the 
iteration mapping 3 has a unique fixed point (^*,\l/*) in /Cm- If necessary, we can further reduce 
o"i > depending only on the data so that wherever a G (0, a\], any fixed point of 3 in /Cm satisfies 
(1237) . □ 



Appendix A. Nonconstant background density 6 for one-dimensional flow 

Proposition A.l. Fix a constant L > 0, and a smooth positive function b(x n ) on [0,L]. Then 
there exists a set of parameters C R 2 x R + so that for any ($ en , £$o,p ex ) £ ^3q ; i/ie system of 
(|1.12p and (|1.13p m TV wit/i f/ie boundary conditions (|1.17p - (|1.2ip /ias a unique one- dimensional 
smooth solution (<p, $) satisfying |V<^| 2 ) > p > and (jl.39p on [0,L] /or positive constants p 
and v\ depending on & en , &o,p ex , b and L. 

Proof. Fix a constant Jo > 0, and we first consider the initial value problem (|1.33j) on the interval 
(0,L) with bo being replaced by a nonconstant positive function b. The initial data (pq,Eq) will be 
specified later. If H'{p) ^ on (0,L) for H'{p) in (fP2l) . (fL33l) is equivalent to 

(A.l) r pVCpWS' (p,E)(0) = (p ,E o ), u = -°. 

[E' = p-b, P 

Claim: For any given constant Jq > 0, there exists a nonempty set ^\{Jq) C R + x R so that for 
any (po,E ) G ^3* (Jo), the initial value problem (jA.lj) has a unique C 1 solution (p,E)(x n ) on the 
interval [0, L] with satisfying p > and the inequality (ll.34p on [0, L] for some positive constant 
v\. Here, the choice of v\ depends on po,E ,J and L. 

Once this claim is proved, one can repeat the argument in the proof of Proposition ll.2l to complete 
the proof of Proposition IA.11 So we only need to prove the claim. 

Fix Jo > and po > supp ^j b + eo satisfying p'(po) — (j^) 2 > v\ > for positive constants eo 
and v\. We will show that if Eq > e± for some e\ > 0, then the initial value problem (jA.ip has a 
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unique G 1 solution (p,E) on the interval [0,L]. Fix a constant Eq > e\. Since p'(p) is C , (jA.ip 
has a unique G 1 solution (p, £7) on [0, e] for a constant e > so that the inequalities 

p>supo+^, £ > ^ and p'(p ) - ( — ) >^ 
[o,L] 2 2 VPo/ 2 

hold on [0,e] where e depends on Eq. Using two equations in (jA.ip . one can directly check that 

, . JL, (pp"(i>) + 2(ff)E 

So if p,E and p'(p) — (^) 2 satisfy the three inequalities above on the [0,e], then we get //(e) > 0, 
E'(e) > and (p'(p) - (y) 2 )'(e) > so (jA.ip is uniquely solvable beyond the interval [0, e], 
moreover (p, E) are increasing functions. Therefore, we can uniquely solve (|A,1|) up to x n = L, and 
the solution (p,E) satisfy p > po, E > Eq and p'(p) — (^) 2 > v\ on [0,L]. For each Jo > 0, we 

have shown that there exists a nonempty parameter set ^I(Jq) so that, for any (pq, Eq) G *P*(Jo), 

(jA.ip has a unique solution (p, E) on [0, J] with satisfying p > and (|1 .39|) on [0, L]. Since ^\{Jq) 

is nonempty, we can find the maximal set ^\{Jq) such that, for any (pq,Eq) G *PJ(Jo), (jA.ip has 

a unique C 1 solution on the interval [0, L] with satisfying p > and (jl.39p on [0,L]. Set ^ = 

U Jo). We repeat the argument in Proposition 1 1 . 21 to get a set ^q of parameters ($ en , ^cbPeai) 
Jo>0 

by using a one-to-one and onto correspondence between (pq,Eq) G and (<3? en , 3§Q,p ex ). This 
completes the proof. □ 

Appendix B. Perturbation of nozzle 

For n > 2, let Ao be an open, connected and bounded subset of R n_1 with G 3 boundary 8Aq, 
and set Mq to be an n-dimensional nozzle defined by 

(B.l) Mq := A x (0, L). 

For < a < 1, let G : A x (-L,2L) R ra_1 be a G 2 ' Q mapping such that G(-,x n ) : Ao — > 
G(Ao x {x n }) C M n_1 is a diffeomorphism for every x n G (— L,2L) satisfying 

(B.2) ||G(-,x n )-/d n _i|| 2 , a , Ao <a 

for a small constant a < which will be specified later. For such a mapping G, we define a 
perturbation of A/"o by 

(B.3) M G := {(G(x),x n ) E R" : x = (x',x n ),x' G A ,x n G (0,L)}. 

A/g is a small perturbation of Mq in the sense that the wall boundary G((9Ao x (0, J)) of Mg 
is a perturbation of the wall boundary <9Ao x (0, J) of Mq- Let T denote the transformation 
T(x',x n ) = (G(x',x n ),x n ) =: (j',y n )- 

Theorem 2. Lei Mq be defined by (jB.ip . Fix 6o > 0, J > and any a G (0,1), and /ei i/ie 
parameter set ^q be as in Proposition Given (Q e n,o, ^o,o,Pex,o) G ^Po> let (^o^o) be the 

corresponding background solution. We assume that 6, (& e n,&ex,Pex) an d are given small per- 
turbations of ($en,0) ®,Pex,o) £ ^Po a nd 3§Qfi, respectively, in the following sense 

\\b - bo\\a,M G < u i 

(B.4) \\§en - ^en,o||2,a,T(r ) + ll^ea; h,a,T(T L ) + \\Pex ~ Pex,0 \\a,T(T L ) < a > 

\&o - &o,o\ ^ *- 
Also, we assume that <l? en and <& ex satisfy 

(B.5) d nw $ en = on T(Tq~ D IQ, 9 ni „l> ea; = on T(r7nT\:) 
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for inward unit normal vector n w on T(8Aq x (— L, 2L)), and that the perturbation T satisfies (|B.2|) 
and 

(b.6) D x T = i n on r ur L . 

Then, there exists a constant a > depending on the data and a such that wherever a € (0, a], the 
nonlinear boundary value problem 

, nr > (div(p(<S>,\V<p\ 2 )Vip) = . 

B.7 < , , in A/g 

1 ' |A<l> = p(<I>,|V ¥ >| 2 )-& 

(B.8) ^ = 0on T , (V^ 2 )) = Pex on T(T L ) 

(B9) * = + onT(r ) 

+ on T(T L ) 

(B.10) = a nro $ = on T(T W ). 

has a unique solution (</?,<&) £ [C^_ a 7-(r )ur(r i ))(-^' G ')] 2 satisfying the following properties: 

(a) The equations in (|B.7|) /orm a uniformly elliptic system in Nq. Equivalently, the solution 
( f, <&) satisfies the inequality 

p\ P ^,\V<p\ 2 ))-\V<p\ 2 >P>0 zn A/^ 

/or a positive constant v, i.e., the flow governed by (<p,$) is subsonic; 

(b) (<p, <&) satisfy the estimate 

m-m II,/, n(-i-a,rcr )ur(r£)) „ ffi , ||(-i-a,T(r<,)ur(ri,)) <r 

/or a m (|B.4|) . The constants u and C depend only on bo, L, $ enj o, &ofi,Pex,o, n, A and a. 
Suppose that (</?, <3>) is a solutions to (|B.7j) - (lB.10p . and let us set 

0(x) :=^oT(x), W(x) := $ o T(x) for x G A/" . 
Then straightforward computation shows that (</>, W) satisfy 
[div x ^4i(x, W, V0, J r ) = 



(B.12) 



\diy x A 2 ^,VW,J T ) = ^(p(W,\J T VcP\ 2 )-b) A " 
.Ai(x,W,V0, Jr) =^2(x,VW, J r ) -nw = on T w 
p(p(W,\J T V ( j ) \ 2 ))=p ex on T L 



for J7-, and A2 defined by 



n 

y=r(x)/i,j = l' 



T 2 M T M T t . ,M T Mq^T 



^(x,z, qi ,M) := ( p (,JJlf q Tf ) __ q T)' i A 2 (x )q2 ,M) := (-^jf-) 
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for (x, z, qi, q2, M) G Wo x R x R n x R n x GL n . Here div x denotes (d Xl , • • • , d Xn )-, and qi and q 2 
denote row vectors. And, n w is inward unit normal vector on T w in (|B.12|) . We rewrite (|B.12p as 

(B.13) 

'div x .4l(x, W, V0, I n ) = div x (A(x, W, V</>, In) - ^l(x, W, Vcj), J T )) 



(=:Hi) . ., 

div x ,4 2 (x, VW, I n ) = pm £Zr 2) ~ b + diVx Wx, VW, I.) - A 2 (x, VW, J T ) ) ° 

("Ha) 

f^i(x,W,V0,4)-n to = (^ 1 (x,W,V^J n )-^li(x,W,V0,J r ))-n to =:s 1 -n to 
\^ 2 (x,V^,J r )-n w = (^ 2 (x,V^/ n )-^ 2 (x,VW r ,Jr))-n u; =:g 2 -n u , ° n " 

p(p(W, |V^| 2 )) = + (p(p(W, |V0| 2 )) - p(p(W, |J r V0| 2 ))) =: p ex + 53 on T L 



so that the left-hand sides of equations in (1B.13P are same as the left-hand sides of (11.24p . and the 
boundary conditions in (lB.13j) can be rewritten as conormal boundary conditions. Also, if ((f), W) 
are sufficiently close to (yo>^o) i n ^(-i~ a T ur L )^ ^ then, by (jB.2|) . there exists a constant C 
depending only on the data and a so that Hi,H 2 ,gi,g 2 and g3 satisfy 

2 3 

(b.14) y: \mi7£ UTL) + E n^iiS;i ouri) < c \% - jtWm < ca. 

1=1 k=l 

The condition ()B.6P implies that Jj- = I n on To U so we have Hi = = H 2 on To U T^, then 
Hi and H 2 satisfy the compatibility condition (|3.ip on To U T^. Now we can repeat the argument 
in £}3]and 2] to prove Theorem [2j 
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